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Abstract 

For a stopped diffusion process in a multidimensional time-dependent domain V, 
we propose and analyse a new procedure consisting in simulating the process with 
an Euler scheme with step size A and stopping it at discrete times (iA)jgN* in a 
modified domain, whose boundary has been appropriately shifted. The shift is locally 
in the direction of the inward normal n(t, x) at any point (t, x) on the parabolic 
boundary of D, and its amplitude is equal to 0.5826(...)|re*o"|(t, x)^f~K where a stands 
for the diffusion coefficient of the process. The procedure is thus extremely easy to 
use. In addition, we prove that the rate of convergence w.r.t. A for the associated 
weak error is higher than without shifting, generalizing previous results by [BGK97] 
obtained for the one dimensional Brownian motion. For this, we establish in full 
generality the asymptotics of the triplet exit time/exit position/overshoot for the 
discretely stopped Euler scheme. Here, the overshoot means the distance to the 
boundary of the process when it exits the domain. Numerical experiments support 
these results. 
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1 Introduction 



1.1 Statement of the problem 

We consider a <i-dimensional diffusion process whose dynamics is given by 

X t = x + ['bis, X s )ds + X s )dW s (1.1) 

JO JO 

where IT is a standard rf'-dimensional Brownian motion defined on a filtered 
probability space (^t)t>0)P) satisfying the usual conditions. The map- 

pings b and a are Lipschitz continuous in space and locally bounded in time, so 
that f ll.ip has a unique strong solution. We consider (_D t ) t > , a time-dependent 
family of smooth bounded domains of R d , that is also smooth with respect 
to t (we refer to paragraph II .5.21 for a precise definition). See Figure [TJ For a 
fixed deterministic time T > 0, this defines a time-space domain 

V= (J {t} x D t = {(*, x) : < t < T, x 6 D t } c]0, T[xE d . 

0<t<T 

Cylindrical domains are specific cases of time-dependent domains of the form 
T> —]0,T[xD, where D is a usual domain of M d (D t = D for any t). Time- 
dependent domains in dimension d = 1 are typically of the form T> — {(t,x) : 
< t < T, <fi(t) < x < (p2(t)} for two functions <p\ and <p 2 (the time- varying 
boundaries). 




Figure 1. Time space domain and its time-sections. 

Now, set t := inf{t > : X t £ D t }, then r A T is the first exit time of (s, X s ) s 
from the time-space domain V. Given continuous functions g, f, k : T> — > R, 
we are interested in estimating the quantity 

ptAT rs 

E x [g( T AT,X TAT )Z TAT + / ZJ(s,X s )ds],Z s = exp(- k(r,X r )dr), (1.2) 

jo jo 

where as usual K x [.] := E[.|X = x] (resp. P z [.] := P[.|X = x]). The approxi- 
mation of such quantities is a well known issue in finance, since it represents in 
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this framework the price of a barrier option, see e.g. Andersen and Brotherton- 
Ratcliffe |ABR96| . These quantities also arise through the Feynman-Kac rep- 
resentation of the solution of a parabolic PDE with Cauchy-Dirichlet boundary 
conditions, see Costantini et al. |CGK06| . They can therefore also be related 
to problems of heat diffusion in time-dependent domains. 

We then choose to approximate the expectation in (11.21) by Monte Carlo sim- 
ulation. This approach is natural and especially relevant compared to deter- 
ministic methods if the dimension d is large. To this end we approximate the 
diffusion ( 11. lj ) by its Euler scheme with time step A > and discretization 
times (ti = iA = iT/m)i> (m G N* so that t m = T). For t > 0, define 
4>(t) = U for ti < t < t i+ i and introduce 

X A =x + j\cj ) (s),X A s) )ds + j\(<j>(s),X$ [s) )dW s . (1.3) 

We now associate to ( 11 M the discrete exit time r A := inf{tj > : X A £ D ti }. 
Approximating the functional V T := g(r A T, X tAT )Z tAT + J rAT Z s f(s, X s ) ds 
by 

rr A AT 

V£ := g(r A A T, A r A AAT )Z r A A/xT + jf Z^f^s), X* s) )ds 
with Zt = e-^ r) ^ )d \ 
we introduce the quantity 

Err (T, A, g, f, k, x) = E X [V T A * - V T ] (1.4) 
that will be referred to as the weak error. 

Note that in V^a, on {r A < T} g is a.s. not evaluated on the side part 
Uo<t<-rW x 9D t of the boundary (g must be understood as a function defined 
in a neighborhood of the boundary). At first sight, this approximation can 
seem coarse. Anyhow, it does not affect the convergence rate and really reduces 
the computational cost with respect to the alternative that would consist in 
taking the projection on &D. It is a commonly observed phenomenon that 
the error is positive when g is positive (overestimation of E X (V T )), because 
we neglect the possible exits between two discrete times: see Boyle and Lau 
BL94| , Baldi |Bal95| . Gobet and Menozzi |GM04| . In addition, it is known 



that the error is of order A 1 / 2 : see |GM04| for lower bound results, see |GM07 



for upper bounds in the more general case of ltd processes. But so far, the 
derivation of an error expansion E^V^a — V T ] = C\fK + o(v / A) had not been 
established: this is one of the intermediary results of the current work (see 
Theorem [4]). 

Our goal goes beyond this result, by designing a simple and very efficient 
improved procedure. We propose to stop the Euler scheme at its exit of a 
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Figure 2. The boundary dD t and the smaller domain D^. 

smaller domain in order to compensate the underestimation of exits and to 
achieve an error of order o(y/A). The smaller domain is defined by its time- 
section 

D A = {xeD t : d{x,dD t ) > coVA\n*a{t,x)\} 

where n(t, x) is the inward normal vector at the closest point of x on the 
boundary dD t , see Figures [2] and [3] for detailfl We shall interpret \n*a(t, x) | as 
the noise amplitude along the normal direction to the boundary. The constant 
cq is defined later in (12.1ft and equals approximatively 0.5826(. . .). Thus, the 
associated exit time of the Euler scheme is given by 

f A = inf{t, > : X* * At) < r A . 

The new Monte Carlo scheme consists in simulating independent realizations 
of 

rr A AT 

V£ = g(f A A T, X*k AT )Z? AAT + J o Z£ is) f(<f>(s),x£ {s) )ds 

and averaging them out to get an estimator of the required quantity E X (V T ). 
Our main result (Theorem [5} is that the asymptotic bias w.r.t. A is signifi- 
cantly improved: 

E x \y&-v T ) = o(VZ) 

(instead of C\f~K + o(y/A) before). This improvement has been already es- 
tablished in the case of the one-dimensional Brownian motion |BGK97| in the 
context of computational finance, exploiting heavily the connection with Gaus- 
sian random walks and some explicit computations available in the Brownian 
motion case. 



1 the closest point of x may not be unique for points x far from dD t . But since 
the above definition of involves only points close to the boundary, this does not 
make any difference. 
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1.2 Contribution of the paper 

To achieve the results in the current very general framework, we combine 
several ingredients (which correspond to the main steps of the proofs). 

(1) We first expand the error E x [V r A — V T ] related to the use of the dis- 
crete Euler scheme in the domain V. Although this issue deserved many 
studies in the literature, the expansion results are new. We prove that it 
relies on the study of the weak convergence of the triplet (exit time, posi- 
tion at exit time, renormalized overshoot at exit time), that is (t a ,X^ a , 
A~ 1 / 2 d(X^ A , 8D t a)) 1 as A goes to 0. This weak convergence result is cru- 
cial in this work and it is new (see Theorem [3j). 

Then, combining this with sharp techniques of error analysis, we derive 
an expansion of the form Err(T, A, g, f, k, x) = C\^K + o(A) in the very 
general framework of stopped diffusions in time-dependent domains. 

(2) Second, we analyse the impact of the boundary shifting, in the continuous 
time problem (see paragraph 12.3.21) . This is related to the differentiability 
ofE x (V T ) w.r.t. the boundary and it has been addressed in |CGK06| . We 
apply directly their results. Then, we obtain the global error estimate of 
the boundary correction procedure (Theorem [5]) . 

We mention that the previous results about the error expansion and correction 
still hold in the stationary setting, see Section [U which also seems to be new. 
A numerical application is discussed in Section [5l Complementary tests are 
presented in |Gob09| . showing that the boundary correction procedure is very 
generic and seems to work without Markovian property for X . This feature 
will be investigated in further research. 

Let us finally mention that we could also consider the diffusion process dis- 
cretely stopped: expansion and correction results below would remain the 
same. 

1.3 Comparison with results in literature 

Up to now, the behavior of (11 A\ had mainly been analysed for cylindrical 
domains, in the killed case, without source and potential terms (i.e. when 
the error writes Err(T, A, g, 0, 0, x) = E x {g(X£)l T & >T ] - E x [g(X T )l T>T ]). Let 
us first mention the work of Broadie et al. |BGK97| , who first derived the 
boundary shifting procedure in the one dimensional geometric Brownian mo- 
tion setting (Black and Scholes model). In |Gob00| and |GM04| . it had been 
shown that, under some (hypo)ellipticity conditions on the coefficients and 
some smoothness of the domain and the coefficients, Err(T, A, g, 0, 0, x) was 
lower and upper bounded at order 1/2 w.r.t. the time-step A. Also, an expan- 
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sion result for the killed Brownian motion in a cone as well as the associated 
correction procedure are available in |Men06j. 

All these works emphasize that the crucial quantity to analyse in order to 
obtain an expansion is the overshoot above the spatial boundary of the dis- 
crete process. In the Brownian one-dimensional framework such analysis goes 
back to Siegmund |Sie79| and Siegmund and Yuh [SY82J. Also a non linear 
renewal theory for random walk, i.e. for a curved boundary, had been de- 
veloped by Siegmund and al, see |Sie85| and references therein, Woodroofe 
|Woo82| and Zhang [Zha88j. We manage to extend their results to obtain the 
asymptotic distribution of the overshoot of the Euler scheme, see Sections [2] 
andGS Concerning the asymptotics of the overshoot of stochastic processes, let 
us mention the works of Alsmeyer |Als94| or Fuh and Lai |FL01| for ergodic 
Markov chains and Doney and Kyprianou for Levy processes |DK06| . These 
works are all based on renewal arguments. 

Finally, for simulating stopped diffusions we also mention the alternative tech- 
nique based on Random Walks on Spheres. This method allows to derive a 
bound for the weak error associated to the approximation of K[V T ] in the ellip- 
tic setting for a cylindrical domain, see Milstein |Mil97| . The same approach 
has also been exploited to obtain some strong error or pathwise bounds for a 
bounded time-space cylindrical domain, see Milstein and Tretyakov |MT99| . 
Recently, Deaconu and Lejay |DL06| have developed similar algorithms, but 
based on random walks on rectangles. However, computationally speaking, our 
approach is presumably more direct. 

1.4 Outline of the paper 

Notations and assumptions used throughout the paper are stated in Section 
11.51 In Section [2] we give our main results concerning the asymptotics of the 
overshoot, the error expansion and the boundary correction. These results 
are proved in Section [3j which is the technical core of the paper. Eventually, 
Section [4] deals with the stationary extension of our results. We still manage to 
obtain an expansion and a correction for elliptic PDEs. Some technical results 
are postponed to the Appendix. 

1.5 General notation and assumptions 
1.5.1 Miscellaneous 

• Differentiation. For smooth functions g(t,x), we denote by d@g(t,x) the 
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derivative of g w.r.t. x according to the multi-index (3, whereas the time- 
derivative of g is denoted by d t g(t,x). The notation Wg(t,x) stands for the 
usual gradient w.r.t. x (as a row vector) and the Hessian matrix of g (w.r.t. 
the space variable x) is denoted by Hg(t,x). 

The second order linear operator L t below stands for the infinitesimal gener- 
ator of the diffusion process X in (11.11) at time t : 

L t g(t,x) = Vg{t,x)b{t,x) + ±Tr(Hg{t,x)[<ra*](t,x)). (1.5) 



• Metric. The Euclidean norm is denoted by | • |. 

We set Bd(x, e) for the usual Euclidean d-dimensional open ball with center x 
and radius e and d(x, C) for the Euclidean distance of a point x to a closed set 
C. The r-neighborhood of C is denoted by Vc(r) = {x : d(x, C) < r} (r > 0). 

• Functions. For an open set V C R x R d and I G N, C^' l (T>') (resp. 
CLiJ'^P 7 )) is the space of continuous functions / defined on V with con- 
tinuous derivatives d^cP t f for \/3\ + 2j < I (resp. defined in a neighborhood 
of V>). Also, for a = I + 9, 9 G]0, 1], / E N, we denote by H a (V) (resp. 
H a (T>')) the Banach space of functions ofC^ l (V') (resp. C^' l {V) ) having 
/ th space derivatives uniformly ^-Holder continuous and |_Z/2j time-derivatives 
uniformly (a/2 — [//2j)-H61der continuous, see Lieberman |Lie96| . p. 46 for 
details. We may simply write C^' 1 or H a when V = R x R d . 

• Floating constants. As usual, we use the same symbol C for all finite, non- 
negative constants which appear in our computations : they may depend on 
V, T, 6, a, g, f, k but they will not depend on A or x. We reserve the notation 
c for constants also independent of T, g, f and k. Other possible dependences 
will be explicitly indicated. 

In the following O po i(A) (resp. 0(A)) stands for every quantity R(A) such 
that, for any k G N one has |-R(A)| < C k A k (resp. |-R(A)| < CA) for a 
constant Ck > (uniformly in the starting point x). 



1.5.2 Time-space domains 

Below, we introduce some usual notations for such domains (see e.g. |Fri64| . 
|Lie96| ). In what follows, for any t > 0, D t is a non empty bounded domain 
of R d , that coincides with the interior of its closure (see |Fri64| , Section 3.2). 
We then define the time-space domain by V := Uo<«tW x A c]0,T[xM d , 
see Figure [U 

Regularity assumptions on the domain V will be formulated in terms of Holder 
spaces with time-space variables (see |Lie96| p. 46 and |Fri64| Section 3.2). 
Namely, we say that the domain V is of class H a , a > 1 if for every boundary 
point (t ,x ) G Uo<*<tO} x 9D t , there exists a neighborhood }t — £g,io + 
el[xBd(xo,eo), an index 1 < i < d and a function <p G H a (]t — ^o^o + 
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Eq[xBiI-i{{Xq, X 1 q , Xq +1 , Xq), £q) S.t. 

{u < t <T {*} x dD t }n{]t -el,t + el[xB d (x ,£o)} 
:= {(*, x) G (]t - 4 t + eg[n[0, T]) x £ d (>o, e ) : 

Xj = ip (t, Xx, Xrf)}. 

If "D is of class H 2 , all domains D t , for t G [0, T], satisfy the uniform interior 
and exterior sphere condition with the same radius r > 0. Moreover, the 
signed spatial distance F, given by 



F{t,x) 



-d(x, dD t ), for x G D c t , d(x, dD t ) < r , < t < T, 
d(x, dD t ), for x G D t , d(x, dD t ) < r , < t < T, 



belongs to H 2 ({(t,x) : < t < T, d(x,dD t ) < r }) (see |Lie96j . Section X.3) 
and n(t,x) = [VF]*(t,x) is the unit inward normal vector to D t at 7ig Dt (x) 
the nearest point to x in dD t (see Figure [3]). The function F can be extended 
as a H 2 ([0,T] x M d ) function, preserving the sign (see |Lie96| . Section X.3). 



1.5.3 Diffusion processes stopped at the boundary 

We specify the properties of the coefficients (b, a) in ( 11.11) with assumption 

(A fl ) (with 9 G]0, 1]) 

1. Smoothness. The functions b and a are in H 1+ g. 

2. Uniform ellipticity. For some ao > 0, it holds £*[crcr*](t, x)£ > ao|^| 2 for 
any (t,x,g) G [0,T] x M d x M d . 

We mention that the additional smoothness of b and a w.r.t. the time variable 
is required for the connection with PDEs. We also introduce assumption (A ) 
for which 2. is replaced by the weaker assumption 

2'. Uniform non characteristic boundary. For some r > there exists ao > 
s.t. VF(t,x)[ao-*)(t,x)VF(t,x)* > a for any (t,x) G Uo<t<T{*} x 
V 9 D t (r ). 

The asymptotic results concerning the overshoot hold true under (Ag), see 
Section 12 . 1L In the following we use the superscript t, x to indicate the usual 
Markovian dependence, i.e. Vs > t, X l ; x = x+ f t s b(u, X^ x )du+ f t s cr(u, X^ x )dW u . 
Now let 

r l ' x := inf{s > t : X 1 / <£ D s } (1.6) 
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be the first exit time of X l s ,x from D s . For functional of the process X stopped 
at the exit from T>, of the form 



u 



(t, x) =E [g(r^ x A T, Xj&^e" fc ( r ' x '")< fr 

+ y e~Jt k ^*r )*/( s> jf*' x )ds], 



(1.7) 



we now recall (see |CGK06| ) that the Feynman-Kac representation holds in the 
time-space domain. Introduce the parabolic boundary VT> = dV\[{0} x Do]. 

Proposition 1 [Feynman-Kac's formula and a priori estimates on u] 



Assume (A*), D £ H b J: 6 H 9j / G H e and # G C°-° with 9 G]0, 1[. Tnen, 
i/iere is a unique solution in C 1,2 (T>) r\C°'°(V) to 



and it is given by (11.71) . 

In addition, if for some 6 g]0, l[,T> is of class Hi+e, 9 G H i+ e t/ien u G Hi + g. 
in particular Vn exists and is 6 -Holder continuous up to the boundary. 
Eventually, for T> G H 3+ e, fc, / G H 1+ g, g G H 3+ e satisfying the first order 
compatibility condition (d t + — k)g(T,x) + /(T, x)| xe gD T = 0, then the 
function u belongs to H 3+ g. 

Proof. The first two existence and uniqueness result for (11.8P are respectively 
implied by Theorems 5.9 and 5.10 and Theorem 6.45 in Lieberman, |Lie96| . 
The probabilistic representation is then a usual verification argument, see e.g. 
Appendix B.l in |CGK06| . The additional smoothness can be derived from 
exercise 4.5 Chapter IV in |Lie96| or Theorem 12, Chapter 3 in |Fri64| . □ 



2 Main Results 

2.1 Controls concerning the overshoot 

The overshoot is the distance of the discretely killed process to the boundary, 
when it exits the domain by its side. To be precise, we use F the signed 
distance function and we consider the quantity F(ij, X A ). It remains positive 
for ti < r A , and at time tj = r A , it becomes non positive. Additionally, 
under the ellipticity assumption, the above inequality is strict: F(t a , X a a ) < 
a.s.. The overshoot is thus defined by F~(r A ,X A A ). Also, since F is in H 2 



d t u + L t u — ku + / = in V, 



(1.8) 



< 



u = g 



on W, 
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(and therefore Lipschitz continuous in time and space), it is easy to see that 
F~(r A , X^a) is of order y/A (in L p -norm for instance). Thus, it is natural to 
study the asymptotics of the rescaled overshoot 



Adapting the proof of Proposition 6 in |GM04| to our time-dependent context, 
see also the proof of Proposition [15] for a simpler version, one has the following 
proposition. 

Proposition 2 (Tightness of the overshoot) Assume (A' e ), and that V 
is of class H 2 . Then, for some c > one has 

sup E x [exp(c[A~ 1/2 F^(s A t a ,X a t a)] 2 )] < +oo. 

A>0,sG[0,T] 



It is quite plain to prove, by pathwise convergence of X A towards X on com- 
pact sets, that (r A AT, X a AaT ) converges in probability to (r AT, X rAT ). The 
next theorem also includes the rescaled overshoot. 

Theorem 3 (Joint limit laws associated to the overshoot) Assume (A' g ). 
and that V is of class H 2 . Let ip be a continuous function with compact support. 
For allt E [0,T], x G Dq, y>0, 



E x [l r A< t Z A Ay2(A^)l F _ (rA>X A A) > yv ^] 



E 



l T < t ZMX T )(l - H(y/\VFa(r,X T )\) 



with H(y) := (Eq^])" 1 /„" P [s r + > z]dz and s := 0,Vn > l,s n := £™ =1 G\ 
the G l being i.i.d. standard centered normal variables, t + := ini{n > : s n > 
0}. 

In other words, (r A , X^ A , A~ 1 ^ 2 F~(t a , X^ a )) weakly converges to 
(t,X t , I VFir(r, X t )\Y) where Y is a random variable independent of (r,X T ), 
and which cumulative function is equal to H. Actually, Y has the asymptotic 
law of the renormalized Brownian overshoot. In the following analysis, the 
mean of the overshoot is an important quantity and it is worth noting that 

one has E(F) = 2E ° ^J""^] := Cq. One knows from |Sie79| that 

C(l/2) , , 

c = - ' = 0.5826... 2.1 

The above theorem is the crucial tool in the derivation of our main results. 
The proof is given in Section I3.ll 
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2.2 Error expansion and boundary correction 



For notational convenience introduce for x £ D { 



'o, 



rrAT 



u(V) = E x (g(r A T, X ThT )Z ThT + 



/ ZJ(s,X s )ds), 
Jo 



/(D) = E I (g(r A AT,I T A , 




Theorem 4 (First order expansion) Under (Ag), for a domain of class 
H2, g £ Hi + 6i, k, f £ Hi+e and for A small enough 

Err (T, A, g, f,k,x) = u A (V) - u(V) 

= c VAE x (l T < T Z T (\7u - Wg)(r,X T ) ■ VF(r,X T )\VFa(r,X T )\) + o(v / A), 
where Cq is defined in fl2.ll) . 

Define now a smaller domain D A C V, which time-section is given by Df = 
{x G D t : d(x,dD t ) > c v^A| VF<r(t, x)|}, see Figure [2l Introduce the exit 
time of the Euler scheme from this smaller domain: f A = inf{tj > : X A G" 
-D A } < r A . The boundary correction procedure consists in simulating 



As above, we do not compute any projection on the boundary. We denote the 
expectation of (JZ2D by u A (V A ). One has: 

Theorem 5 (Boundary correction) Under the assumptions of Theorem^ 
if we additionally suppose V-F(., .)|V-Fa(., .)| is in C ,2 ; then one has: 



The additional assumption is due to technical considerations to ensure that 
the modified domain V A is also of class H 2 . It is automatically fulfilled for 
domains of class C 3 and a in C 1 ' 2 . 

2.3 Proof of Theorems^ and [5] 

2.3.1 Error expansion 

By usual weak convergence arguments, Theorem 0] is a direct consequence 
of Proposition [2] (tightness), Theorem [3] (joint limit laws associated to the 




(2.2) 



u A (V A ) -u{V) 



o(VA). 
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overshoot) and Theorem [6] below. 



Theorem 6 (First order approximation) Under the assumptions of The- 
orem^ one has 

u A {V) - u{V) = o(v / A) + 

E x (l T ^ T Z A A (Vu - Vg)(r A ,iT dDTA (X^)) ■ VF(r A , X A A )F"(r A , X A A )). 

Remark 7 In the above statement, we use projections on a non convex set, 
which needs a clarification. With the notation of Section \1.5.2. introduce T r ° : = 
inf{s > : X A ^ Vd s (Vo)}- For s G [0, r r °] the projection 7i£) s (X A ) is uniquely 
defined by 



*D S (X A ) = X? + (VF)*(,,X A )F-( S ,X A ), 



(2.3) 



see Figure^ Large deviation arguments (see Lemma\Bbelow) also give F x [r r ° < 
t a < T] = Op i(A). Thus, in the following, for s > r r ° , 7rg s (X A ) and 
ttq£, s (X a ) denote an arbitrary point on dD s . This choice yields an exponen- 
tially small contribution in our estimates. 



x 



>4 



F-(t,x) 



^dD t (x) = y 



\ n{t,y) = [VF]*(t,x) 



Figure 3. Orthogonal projection ^dD t { x ) of x ^ D t onto the bound- 
ary dDt and the related signed distance F(t,x). Here F(t,x) < and 
d(x,dD t ) = \F(t,x)\=F-(t,x). 



Proof. Denote e A := u A (V) — u(V) the above error. Write now 

e A =E x [g(r A A T, X A A/sT )Z A A/sT - g(r A A T, 7r fl ^ Ar (X£ AT ))Z£ AT ] 

rr A AT 

+ {E x [g(r A A T, tt^JX^t^zXt + / Z A s) f(<p(s), X A s) )ds] 

J 

-u(0,X A )} 
:=e A + e A 
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We introduce here the projection for the error analysis. From fl 2 . 3 H and Propo- 
sition El a Taylor expansion yields 

ef = - E x [l T ^ T Z^Vg(r A ,n dDTA (X^)) • VF(r A X a a )F~(t a , X A A )] 

+ 0(A( 1+ ^ 2 ). (2.4) 



In the following, we write U = V (resp 17 < V) when the equality between U 
and V holds in mean up to a O po i(A) (resp. E x (£/) < E^V) + O po i(A)). We 
also use the notation U = 0(V) between two random variables U and V if 
for a constant C, one has \U\ < C\V\. Because g(r A A T, ttjj a t (X^ aT )) = 
u(r A A T, 7T£, A (X^ Ar )), we can write a telescopic summation: 



A E 



£ ^ i+1 ,7r s (X A +1 M 



0<ti<r A AT 

- life, 7T Stj (X A ))Z A + Z A f (t J7 X A ) A) 1 T , 0>T A AT 
0<ti<T 

-n(t,,X A )Z A + Z A f(t 4 ,X A )Al) V 0>t a aT 



since for tj < r A , X A G and thus 7r£> t .(X A ) = X A . To proceed, the 
key idea is to introduce on the event {£; < r A }, the partition {F(ti,X A ) G 
(O^A^ 1 " 6 )]} U {F(t u X A ) > 2A^ 1 " £ )} := A 5 . U (A £ u ) c , e > 0. This allows to 
split the cases for which X A is close or not to the boundary dD ti . Lemma [8] 
ensures that (X A ) se [ ti ,t i+1 ] stayed in B(X A , A^ 1_e )) with a probability expo- 
nentially close to one. Then, on (A e .) c , the smoothness of the domain yields 
l (Af ) oP[X A +l G A 1+1 |7 r tl ] = 1 - 0(exp(-cA- £ )), see Proposition [TJ for a 

proof of this claim. On the other hand, on A E ., X A is sufficiently close to the 
boundary to make the contribution of the overshoot at time tj+i significant 
for the error analysis. Write: 



A E 



0<U<T 



£ _ l tl<T A {l^ 7T Bti+i (X A + J)Z A +1 

^,X A )Z A + Z A /(^,X A )A 



+ 1 ( A t l ) - L V S 6^,t I+1 ], XAeB(A t A ,A2 (1 - £) 



M (t m ,X A +i )Z A +i 



-u(t h X A )Z A + Z A f (t,, X A ) A] }) l r , 0>T A AT := e A + e A . (2.5) 



Let us first deal with e A x . In our framework, u is (1 + #)/2-H61der continuous 
in time and V« is ^-Holder continuous in space on a neighborhood of T>. A 
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Taylor expansion at order one and the equality (12.31) give 



4i=( E W^f. z £Vu(t h x£).VF(t i+1 ,x£ +i )F-(t i+1 ,x£ +1 ) 



0<t;<T 



-O(|F-(t t+1 ,X^ i )| 1+9 ) + O(|X4 i -X < t| 1+e ) + O(A^)l)v 0>rAAT 



+ J2 i tj<T Ai^[o(|F-(t i+1 ,x t f +l ; 

o<u<t 

+ 0(|X£ +i - X f t| e F-(t J+1 , X A J) + 0(A 



E ' WZ T A A Vu(/, X A A ) • VF(r A , X A A )F" (r A , X A A ) 

A ^ 1 o(\x£ +i -x£\^ 

l T r a>T A AT 



1 + . 
2 



where we used once again Lemma [8] for the last equality. Standard arguments 
yieldE[|X A i -X A ^|^J = O(Af)foranyp>0andE[|F-(t m ,X A 1 )|P|J- t J = 
E[|F-(t m ,X t f +i ) - F-(U,X£)\v\F t% ] = O(Al) on {U < r A }. Thus, we can 
now rewrite 

e A I(l r A< T Z r A A Vn(r A ,7r aDTA (X r A A )) ■ VF(r A , X A A )F"(r A , X A A )) V 0>T z 



AT 



+ e 211; 



4i=( E W1a?O(A^))v 0>t 



'AT- 



0<ii<T 



To handle e A n the idea is to use the occupation time formula and some sharp 
estimates concerning the local time of (F(s, X^)) s<TAt a in a neighborhood of 
the boundary. We have 



e A n | < CA^[ A 



l+S 



-1 



TAt 



TAt 



LF(<^(t),X^ (t) )e[0,2AV2(l-£)]^ I l r ''0>T A AT 



< CA+ A" 1 

E i_i_n / -' \ 

<CA^ (A- l j_ Aimi s) L y TATA (F(.,X A ))dy V 0>t a aT 



lT(t,X t A )e[-A 1 /2(i-e) j 3A 1 /2(i-e)]^J lr r 0>r A AT 
3A l/2(l- £ ) 



where we have used Lemma[8]at the second equality and the uniform ellipticity 
assumption for the last one. Now an easy adaptation of the proof of Lemma 
17 |GM04| to our time-dependent domain framework gives 

E[^ ArA (F(.,X A ))]<C%| + A^). (2.6) 

jg l-\-6 1 

Thus, one has |e A n | < CA~ 2 = o(Aa) for e small enough. Hence, the above 
estimates and Lemma [8] give 

e A ^ (l T A <r Z A A Vn(T A ,7r^ A (X A A ))-VF(r A ,X A A )^(r A ,X A A ))+o(A5). 

(2.7) 
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Let us now turn to e^. If g G ii-i+e (which implies u G 'H-z+e in view of Propo- 
sition Q]), the term can be handled with somehow standard techniques. 
Namely Taylor like expansions in the spirit of Talay and Tubaro |TT90| . For 
simplicity we handle under the previous smoothness assumption on g and 
u. The proof under weaker assumptions (g G Hi + g), that involves sharp esti- 
mates on possibly exploding derivatives of u near the boundary, is postponed 
to the Appendix. We recall that 



d A E 
"22 ~~ 



0<U<T 



*.<r Ai (4) 



Vse{ti,t i+1 ], x*eB(x t A ,A? {1 ~ E) ) 



u(t i+u x£)z£ +i 



-u(U,x£)z£ + z£f(u,x£)A]) Vo> t a 



For all (s,y) G P introduce the operators L S)y : C 1,2 (P) — > C(V), ip \— > 
i ► L S: yLp(t,x) = Vf(t,x)b(s,y) + ^Ti[HLf(t,x)[aa*](s,y)]). Recalling 
that X£) + L t . jX Au(ti, X%) - ku(U, X£) + f(U, X%) = 0, Ito's formula 

gives 



X A €-B(X t A ,A5 Cl_e) ) 



e| = ( „| : <r 1 " < " 1,Ay0l «'«« 

/'■ + '(Z S A - Z, a )(& + L t x& - k(ti,X£))u(s,X?)ds 

Jti '« 

- (d s + L tuXt - k(t t ,X£))u(t t ,X£)) 



ds 



M, 



L T r l)>T A ATl 



(2.8) 



where for all v G [tj,tj + i], M t v := / Z^Vu(s, X^)a(U, Xp)dW s is a square- 

Jti . . ... 1 

integrable martingale term. Note that in this definition, in whole generality, 

M ti>v is not stopped at the exit time r ti := inf{s > ti : G" -D s }. If r t . < 
(which happens with exponentially small probability on (AlJ°), the term 
Vu(s, ), s G [t^, tj+i] in M tuti+1 has to be understood as the smooth exten- 
sion of V« to the whole space. In particular this extension remains bounded. 
Now, we derive from Lemma [8] 



EU( E 1 * i <r Al (^) 0l V sefe)t . + l]jX A eS(X A iA ^l-e) ) ^ > ii + l 
0<U<T 1 + 1 y '< 



L r r o>T A ATj 



0<ti<T 
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We can thus neglect the contribution of the martingale terms in (12.81) . We now 
develop the other quantities in (12.81) with Taylor integral formulas to derive 



[ ti+ \z* - Z t A )(d s + L ti)X A - k(t u X£))u(s,X*)ds 

£ + \d t u(s,X*) -d t u(ti,X£))ds 

= ( U+ \d t u(t h X£)cr(t u X£)(W s - W u )ds 

+ 0(A 1+1 * £ [9 t n] tii± ,+A 2 |V^n| 0O + A sup |X A - X* \ 1+e [Vd t u] x , e 

2 se[ti,t i+ i] 



rti+i . . 

/ (L t x au(s, X s ) - L t X AU[ti, X t ))ds 

Jti i H 

= f l+ \H u {t u X^)a{t u X^){W s - W u ),b(ti,X£))ds 

+ ^J^ +1 Tv{{DMtuXtM^ X u)(W s - W ti )) ■ a(U,X£))ds 

+ 0(A 2 {| J D 2 u| 00 + \lFu\oo + |9tVu|oc} + A 1+1 -^[D 2 u} tl ±0 

+ A\D 3 u\oo sup |X S A -X A | 2 + A sup \X?-X*\ 1+ *[lPu]J) t 

s£{ti,t i+1 ] s£[ti,t i+1 ] 

k(U,X£) r +1 (n( S ,X A ) - u{t u X£))ds 

J ti 

=k(U,X*) f U+1 Vu(U,X^)a(U,X t %W s - W u )ds 

J ti 

+ 0(A 2 (\d t u\ 00 + \Vu\ 00 ) + A\D 2 u\ 00 sup |X A -X A | 2 ), (2.9) 

s£[ti,ti + i] 



where [•]*,«, [-}x,a, ot £ (0, 1] denote respectively the Holder norms of order a 
in time and space (see Chapter IV Section 1 p. 46 in |Lie96| for a precise 
definition). 



Hence, bringing together our estimates and exploiting the relations between 
the spatial and time derivatives for u (through the PDE) , from ( 12.81) and (12.91) 
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we derive 

A E 



e 22-( W 1 (^.)° 1 V,6h,t (+1 ],XA6B(^,Ai! 1 -)) 

0(A 2 {1 + 



+ 0(A 1+1 ^{1 + \u\oo + (Vuloo + |D 2 u|oo + |i^ 3 u|oo + [D 2 u] tt i±e} 

+ 0(A sup \X A -X A \ 1+9 {l + \u\ 00 + \Vu\ 00 + \D 2 u\ 00 + \D 3 u\ 00 + [D z 
se[ti,t i+1 ] 

+ 0(A sup \X? - XffUlPuloo + | J D 3 u| 00 } 

«6[*i)*<+l] 



+ M ti>ii+1 jjl r , 0>T A AT , (2.10) 

where M ti jii+1 denotes the sum of the terms involving the Brownian increment 
(W s — W ti ) S £{t u t i+1 ] in the above equations (12 .9j) . Under our current assump- 
tion, i.e. it G H 3+0 , all the norms appearing in (12.101) and all the derivatives 
appearing in the (M ti}ti+1 )o< ti< T are bounded. Hence, 

=l ti<T Al {Al) cE[M ti , ti+1 \T ti ] + O pol (A) = O pol (A), (2.11) 



E 



\e£ 2 \ <CA~. (2.12) 

Plug flU]) and ([27121 into (E5|). The statement is derived from fl2j) and (1231) . 
We specify in the Appendix how to complete the proof from a sharper version 
of ( 12TT01) deriving from (JUj), when g G H 1+g . □ 



2.3.2 Boundary Correction 
One has 

u A (D A ) - u(P) = [u A (V A ) - u(V A )] + [u(V A ) - u(V)}. (2.13) 

(1) The first contribution in (12.131) has been previously analysed in Theorem 
[H except that the domain V A depends on A. We can show that it is equal 
to c v / AE(l T < r Z T (Vn - Vg)(r, X T ) • VF(r, X T ) \VFa(r, X T )\) + o(y/A). 
We briefly sketch the proof of this assertion, which is done in two steps. 
For this, set u A = u(V A ) for the solution of the PDE in the domain V A . 
• Step 1. It is well known that all PDE estimates depend only on bounds 
on the derivatives of the level set functions (y2 ) arising in the definition 
of the time-dependent domains (see section fl.5.21) . and on the bounds 
on the derivatives of data g, / and k. Hence, since V A is a small per- 
turbation of class H 2 (because VF\VFa\ has this regularity) of the 
domain V of class H 2 , all PDE estimates on u A remain locally uni- 
form w.r.t. A. In addition, u A and its gradient converge uniformly to 
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u and Vu. This argumentation allows us to state that the first order 
approximation theorem holds: 

u A (V A )-u(V A ) = o(v / A)+ 

E,(^< T 4(V« - Vg)(t*,ir aD *(X&)) ■ VF a (t a ,X a a )[F a ]-(t a ,X a a 

T 

where F A and f A are respectively the signed distance to the side of V A 

and the related discrete exit time. 
• Step 2. The second step is to prove that the analogous version of 

Theorem [3] holds, with f A instead of r A . Actually, a careful reading of 

its proof shows that it is indeed the case, without modification. 
(2) Finally, the last term in (12.131) is related to the sensitivity of a Dirichlet 
problem with respect to the domain. By an application of Theorem 2.2 in 
|CQK06| with Q(t,x) = -c VF(t,x)\VFa(t,x)\ (in C 1 ' 2 ), one gets that 
this contribution equals 

-coVAE(l T <r^r(Vu - Vg)(r,X T ) ■ VF(t,X t )\VF<t(t,X t )\) + o(VK). 
This proves that the new procedure has an error o(y/A). □ 

3 Technical results concerning the overshoot 

This section is devoted to the proof of Theorem [3j We first state some useful 
auxiliary results. 

Lemma 8 (Bernstein's inequality) Assume (Ag-1). Consider two stop- 
ping times S, S' upper bounded by T with 0<S' — S<Q<T. Then for any 
p > 1, there are some constants c > and C := C((Ag-l) , T), such that for 
any r) > 0, one has a.s: 

P[ SUp \X A -X A \ > 7] 

te[s,s'} 

E[ sup |Af -XJT 

te[S,S"] 

For a proof of the first inequality we refer to Chapter 3, §3 in |RY99| . The last 
inequality easily follows from the first one or from the BDG inequalities. 

Lemma 9 (Convergence of exit time) Assume (A' e ) and that the domain 
is of class H 2 . The following convergences hold in probability: 

(1 ) lim A ^o t a AT = t AT; 

( 2) lim A ^o X A a aT = X ThT ; 



Fs] <C exp — c 
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F s ] <ce p/2 . 
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(3) limA-o sup t < T \Xf {t) - X t \ = 0. 

The proof of the first two assertions in the case of space-time domain is analo- 
gous to the case of cylindrical domain (see |GM05| ) and thus left to the reader. 
The last convergence is standard. 

The following results are key tools to prove Theorem [31 A similar version is 
proved in |Sie79| . but here, we additionally prove the uniform convergence. 

Lemma 10 (Asymptotic independence of the overshoot and the dis- 
crete exit time). Let W be a standard one dimensional BM. Put x > and 

consider the domain V :=]0,T[x] — oo,x[. With the notation of Section^ for 
any e > we have 



lim sup P [r A < t, (W t a - x) < yVX] -P [r < t]H(y) 

A >0 t£[Q,T],y>a,x>AV2-z 



= 0. 
(3.1) 



If the Euler scheme starts close to the boundary at a small distance d, its 
discrete exit likely occurs after a time roughly equal to d 2 . This feature is 
quantified in the above lemma. 

Lemma 11 Assume (A' g ). and that the domain is of class H 2 . Let < (3 < 
a < 1/2. For all rj > 0, there exists C := C v > s.t. for A small enough, 
Vs G AN n [0, T] and Vx E V dDs (A a ) n D s , one has 

F[t a A T > A 2/3 |A S A = x}< C{^^ + A"), 

where r A := inf{t; > s : X t A (£ D ti } . 

Lemma 12 Assume (A' e ), and that the domain is of class H 2 . There exists 
C > 0, such that Vs G AN fl [0, T], Vx G D s . Vt G [s, T] and Mb > a > 0, one 
has 

P[r A < t, A" 1 / 2 F^(r A , A A A ) G [a, b] |A S A = x] <C((6 - a) + A 1 / 4 ) 

where r A is shifted as in the previous lemma. 

The proof of these three lemmas is postponed to Section I3.2L 

We mention that if oo* is uniformly elliptic, Lemma [12] is valid without the 
A 1 / 4 (see the proof for details). In that case, it means that the law of the 
renormalized overshoot is absolutely continuous w.r.t. the Lebesgue measure 
on 1R + , with a bounded density. This is also true at the limit, in view of 
Theorem [3j 



19 



3. 1 Proof of Theorem 



Consider first the case V =}0,T[xD where D is a half space. The theorem 
in the case of BM is then a direct consequence of Lemma Now to deal 
with the Euler scheme, we introduce a first neighborhood whose distance to 
the boundary goes to with A at a speed lower than A 1 / 2 (below, the speed 
is tuned by a parameter a, see Figured)). The characteristic exit time for a 
starting point in this neighborhood is short (Lemma [II]) , thus the diffusion 
coefficients are somehow constant and we are almost in the BM framework. 
Also, a second localization w.r.t. to the hitting time of this neighborhood guar- 
antees that up to a rescaling we are far enough from the boundary to apply 
the renewal arguments needed for the asymptotic law of the overshoot (this is 
tuned by another parameter e, see Figured]). 

For a more general time-space domain of class H 2 two additional tools are 
used: a time-space change of chart and a local half space approximation of the 
domain by some tangent hyperplane. 

For notational convenience, we assume from now on that the time-section do- 
mains (-Dt)te[o,T] are convex so that HdD t is always uniquely defined on D^. To 
handle the case of general H 2 domains, an additional localization procedure 
similar to the one of Theorem [6] is needed. We leave it to the reader. 
For the sake of clarity, we also assume k = (Z = 1). This is an easy simpli- 
fication since owing to Lemma [9j Z^. AT converges to Z tA t in L\. 

Step 1: preliminary localization. For a < 1/2 specified later on, define 
r A a := mf{U > : F(U,X£) < A a } < r A . We aim at studying the conver- 
gence of 

& A {t,x,y) : = E x [l T A< tiF - (T A iX A A) > y ^(X A A )] 
and for this, we define for all < s < t < T (s e AN), (x, y) G R d x R+ 

V A (s,t,x,y) :=P[r A < t,F~(r A ,X?k) > yVA\X* = x], 
A(t,a,e) :={r Aa < r A ,r A « < t, F(r Aa , X^J > A 1/2 ~ e }. 

Here, e is a fixed parameter in ]0, l/2[, such that a < 1/2 — e (take e = 
(a + l/2)/2 for instance). 

In the definition of ^a, t A has to be understood as the shifted exit time 
inf{t< > s : X A £ D ti }. By Lemma El ^ x [r A = r A « < t] + F x [r Aa < 
t,F{r A «,X* a ) < A 1/2 ~ e j = O pol (A) using a < 1/2 - e. Hence, 

& A (t,x,y) =E iB [lA(t, a , e ),F-(TA > x^)>»VS¥'( j ^) 1 TA<t] +O po i(A) 

+ E x [l A{t ^ E) ip{X? A a )^ A (r Aa , t, X* a , y)} + O pol (A). 
The first term in the right hand side above converges to 0, using the conver- 
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gence in probability of \X^ AT — X^ aAT \ to (analogously to Lemmal). This 
gives 

V A {t, x, y) = E x [l A{t ^ E)l p(X* J*a(ta«, *, , 2/)] + o(l). (3.2) 




—\ 1 time 

t T 



Figure 4. The two localization neighborhoods with a < — e. 

Let us comment again these two localisations. That with A" enables us to 
freeze the coefficients of the Euler scheme, because the exit time is likely close 
to the initial time. That with A 1 / 2-6 ensures that it starts far enough from the 
boundary to induce the limiting behavior of the overshoot. This right balance 
regarding the distance of the initial point to the boundary is crucial. The final 
choice of a (and thus e) depends on the regularity 6 of the coefficients b and 
a. 

Now, it remains to study the convergence of ^a(-)- 

Step 2: diffusion with frozen coefficients. Denote r Aa := s, X^ a := x. 

Conditionally to introduce now the one dimensional process (Y s ) s >s, Y s = 
F(s,x) + (VFcr)(s, x)(W s — W s ). Note that we do not take into account the 
drift part in the frozen process. From the next localization procedure, it yields 
a negligible term. Since Y has constant coefficients, we apply below LemmafTOl 
to handle the overshoot of Y w.r.t. R + * . Define r A ' y := inf{i; > s : Y ti < 0} 
and rewrite 

y A (s,t,x,y) := V%(s,t,x,y) + R A (s,t,x,y), (3.3) 
H><i{s,t,x,y) := P^[r A ' y < t, (y T A,r)- > yVZ]. 

From (A e -2') that guarantees that Y has a non degenerate variance and 
Lemma [TOJ, one gets 

sup \^1{s,t,x,y) - P M [r A > y < t] ( 1 - H(y / | ( WFa) {§, x ) | ) ) | — > , 

where A a > £ := {(t, x) : < t < T, x G V dDtiA <*)\V aDt (A l/2 ~ E )}- Plug now this 
identity in (I3.3P to obtain with the same uniformity 

* a(s, t, x, y) =P^.[r A ' y < - H(y/\(VFa)(s, x)\)) + R A (s, t, x, y) + o(l). 

(3.4) 

Step 3: control of the rests. We now show that R A (s,t,x,y) = o(l) where 
the rest is still uniform for (s, x) G ^l a ' e . This part is long and technical. First, 
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decomposing the space using the events {r A = r A ' y }, {F (r A , X A A ) > y\/~A}, 
{(Y t a,y)~ > y\/A} and their complementary events, write: 



\R A (s,t,x,y)\ < R A (s,t,x) 

+ P^[r A < t, F-(t a , X a a ) > yVA, (Y T ., Y y < yVX, r A = r A ' y ] 

+ P^[r A < t,F~(r A ,X A A ) < yVE, (Y t a, y )~ > y^,r A = r A ' y ] (3.5) 

with Ri(s,t,x) < P^[r A < t,r A ± r A ' Y ] + P s -,,[r A ' y < t, r A ± r A > y ] : = 
(i?^ 1 + R A )(s,t,x). Let yA be a given positive function of the time-step s.t. 
Va — * specified later on. 

On the event {r A = r A > y , \Y t a,y - F(r A ' y , X A A , y ) | < y A VA}, the conditions 
F-(r A ,X A ) > yv /A and (Y^r)" < yV^" imply A~ 1/2 (Y rA ,y)~ G [j/-3/a,2/). 
Similarly, (Y^r)" > y\/A~ and i^(r A ,X A A ) < y\/A imply A" 1 / 2 ^^)" 
G [y, 7/ + J/a) ■ Hence, by setting 

R 2 A (s,t,x) := 2P^[r A ' y < t, r A = r A ' ¥ , \Y T ., Y - F(r A ' y , X A A , y )| > y A VA], 
R 3 A (sXx,y) :=P^[r A ' y < t, A" 1 / 2 ^^)- G [y - y A , 2/ + 2/a), r A = r A ' y ], 



we obtain R A (s, t, x,y)\ < (R A + R A )(s, t, x) + R A (s, t, x, y). 
Term R A (s,t,x,y). From Lemma [T2l applied to the process with frozen coef- 
ficients, one gets 

R A (s,t,x,y)<C(y A + A 1 / 4 ). (3.6) 

Term R A (s,t,x). Let us explain the leading ideas of the estimates below. 
Usually, it is easy to prove inequalities like \Y t — F(t, X A )\ L2 = 0(A 1 ^ 2 ) (for 
a fixed t), but this not enough to control R\. To achieve our goal, we take 
advantage of the fact that the time t is the stopping time r A,y which is likely 
close to s. Thus, Y t a,y — F(r A,y , X A A}Y ) should be much smaller that A 1//2 in 
L2-norm. 

Introduce for < (3 < a < 1/2, t a p ■= inf{s > s : \X A - x\ > A^} A (s + 
A 5 ), 5 := 2(3 + 7, 7 > 0. Clearly, one has 



T A > Y <t,T A = T A > Y ,T A <T A ,, 



\R 2 A (s : t : x)\ <2P S ~,* 

\Y t a,y - F(r A < y ,X A A , y )| > y A VA\ + 2P s -, £ [r A > r A .,r A < t] 
:=( J R 2 , 1 + J R 22 )(S,t,x). 

Let us first deal with R A (s, t, x). By the Markov inequality, one has 

R%{8,t,x) < 2A- 1 yl 2 E^[l T A <TA ^ T A,Y< t ^ =T A,Y\Y T A, Y - F(r A ' y , X A A , y 



)| 2 ]- 
(3.7) 



Note that since V is of class H 2 , F has the same regularity, i.e. it is uniformly 
Lipschitz continuous in time, its first space derivatives are uniformly Lipschitz 
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continuous in space and l/2-H61der continuous in time. Thus, assuming up to 
a regularization procedure that F G C 1,2 ([0,T] x K d ), Ito's formula yields for 
all t > S, 

F(t, X t A ) =F(S, x) + J!VF{ 8 , X*)dX* 

+ J^d s F(s,X^) + ^Tr(H F ( S ,X^aa*(<f ) (s),X^ ) )))ds 
:=F(§, x) + Nf{s, Xf)a{(t>{s),X^ a) )dW a + R${s, t, x) (3.8) 

J s 

=Y t + R${s, t, x) + f_ (VF(s, X*)a(0(s),X* s) ) - [VFa] (s, x))dW s . 

From (A e -1) and the assumptions on V one derives \R$\(s,t,x) <C{t-s). 
Thus, for any given stopping time U G [s, T A ,s], the working assumptions (i.e. 
smoothness of a, F) and standard computations yield 

E[\F(U,X$) - Y v \ 2 ] < C{A 2 ? +5 + A 8 ^). 

From ( 13 .7R and the above control with U = r A,y A t a p, one obtains 

R%{s,t,x) < Cy A 2 A~\A 2l3+s + A 8 ^). (3.9) 

Let us now control R A (s, t, x). From Lemmas l8l and [TT1 for any rj > we write 

R 2 A 2 (s, t, x) < F s>i [r A0 <S + A S ]+ P s - 5 [r A A t > 5 + A 5 ] 

< C„(exp (-cA 2 ^) + A"""" 5 / 2 + A 5 / 2 ). (3.10) 

Take now a = < 1/2, ^ = 16(g e +1) , 7 = g^y, y A = A 9/16 . Check that 

for 5 = 2(3 + 7 = 2a — 4r/, one has 5 = 1 ^g 4 , /5 = ^W^jy < a, 3r] < a. 

Thus, i?i 2 (5,t,x) = 0(A r '). In addition, y^A 5 ^- 1 = A 6 '\ y^A 2 ^ 8 ' 1 = 
0(A 1 /( 8 ( 1 + 9 ))). Hence, from and rt!UT)D 

t, x) < c(A 1 ^ 1+d ^ + A e / 8 + A^ 16 ^ 1 ") < CA e / 32 . (3.11) 

Term R A (s, t, x). We give an upper bound for R A (s, t, x). The term R A (s, t, x) 
can be handled in the same way. From the previous control on R A (s, t, x) and 
for the previous parameters, one gets 

R A \s, t, x) =P s ,,[r A < t, r A ^ r A ' y , r A < r A „] + 0(A") 
=P^[r A <t,r A >r A ' y ,r A <r A ,] 

+ P^[r A < t, r A < r A > y , r A < r A ,] + 0(A"). 

Then, splitting the first probability according to A~ 1 ^ 2 (Y t a,y)~ < y A or not, 
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and the second one according to A l ' 2 F (t a ,X a a ) < y A or not, we obtain 
R%(s t t,x) 



<(P s -,,[r A ' y <t,A^(Y T A, Y )~ <y A ] 

+ P^[r A < t,r A > r A ' y ,r A < r M , A~^\Y T ., Y - F(t a > y , X A A > y A 
+ (P s -, £ [r A < t,r A < t a > y ,t a < r A „, A^/ 2 |F t a - F(r A ,X A A )| > y A ] 

+ P s -,,[r A < t, A~ 1/2 F~(t a , X%) < y A ])+CA\ 



for the previous function (?/a)a>o- Since we could obtain the same type of 
bound for R A (s,t,x), from Lemma [T2l and following the computations that 
gave (13.91) we derive for the previous set of parameters 

Rtts,t,x) < C{y A 2 A-\A^ +S + A 5 ^) + A r > + y A + A 1 / 4 ) < CA e ' Z2 . 

(3.12) 

From firm (EUl), (ESI) we finally obtain R A (s,t,x,y) = 0(A 9 ^ 2 ) = o(l). 
The rest is uniform w.r.t. (s,x,y) G ^4 a ' e x R + . 

Step 4. Final step. Plug the previous results in (13.41) . We derive from (13.21) 

<P A (t,x,y) = E x [l A{t ^ £) if(X A a ) 

x P ^,^ a [^ A ' y < " ^/|VFa(r AQ ,X A J|))] + o(l). 

Moreover, note that taking y = in the previous controls gives immediately 

P s ~, 5 (r A ' y <t)-P s -, i (r A <t) = o(l) 

uniformly in (s, x) G *4. a ' e . Thus, we finally obtain 

^ A (t,x,y) = E,[l A(t , a , # (I^)l T A< t (l - ff(y/|VFa(r A «,^ a )|))] + o(l). 
Under continuity arguments as in step 1 (localization), we eventually get 

^ A (t,x,y) = E x .[l T A<^(X A A )(l- J f/(y/|VFa(r A ,X A A )|))] + o(l). 
We complete the proof using Lemma M 



Mt,z,y) ^ Q E x [l T <MXr)(l - H(y/\VFa{T,X T )\))}. 



□ 



3. 2 Proof of Lemmas QJjj, Q7] and [7 
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Proof of Lemma [J7J We shall insist on the dependence of the exit times with 
respect to x, by setting r A := inf{tj = iA > : W ti > x} : = t x and 
analogously for r = t x . Our proof relies on the following convergence (see 
equation (19) in Siegmund |Sie79| ): if we set (for any y, z > 0) 



D(z, y) = P [W t a - z < yVA] - H(y), 

then 

lim \D(z,y)\ = 0. 

2 A- 1 / 2 -»+oo 

Using the monotonicity and the uniform continuity of H(y), Dini's Theorem 
yields that the above limit is actually uniform with respect to y > 0. It follows 

sup \D(z,y)\ - 0. (3.13) 

y>0, z£[A 1 i /2 ~ £ / :i ,00) 

Additionnally, we have 

sup |P (r A > t) - P (r x . > t)\ - 0. (3.14) 

x>o, telA 1 - 4 */ 3 ?] 

To prove this, we apply Theorem 3.4 in |Avi07| which states that 
supM\l M<x - 1 A I < 3(sup f M (m)\\M - M\\ Lp )^ 

xGM met 

for any p > and for any random variables M and M, such that M has a 
bounded density /m(-)- Now, consider M = sup s<t W s and M = sup s=iA<4 W s . 
The density of M is bounded by 2/ \j2nt. On the other hand, Lemma 6 in 
|AOP95| gives \\M - M\\ Lp < C P (T)A 1 / 2 . Hence, we get for t > A 1 " 4 ^ 3 , 

|P (r x A > t)-F (r x > t)\ < E\1m <x - W| < C P (T)A^ , 
which leads to (13.141) . 

We can now proceed to the proof of Lemma [TQl assuming that x > A l l 2 ~ £ . 
First, note that if x/\/i > A _e/3 -> +00 as A -> 0, P (r A < t) and P (r x . < t) 
are both O po i(A). Thus, the difference in Lemma [TOl converges to as A — > 0. 



Suppose now that x/y/i < A £ / 3 , hence \/^ > xA £ / 3 > A 1 / 2 2e / 3 , and write 
for t E AN* 

P := P [r A > t, W t a -x< yVA] = / + °V' A (0, x - z)F [W t a - z < yVA]dz 

Jo 

where qt ,A (-, •) denotes the transition density of the Brownian motion dis- 
cretely killed at level x. Introduce the partition R + = [0, A 1/2 ~ £/3 )U[A 1/2_e/3 , +00). 
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Then, 

f + OO 



r+oo 

P = R+ q? A (0,x - z)D(z,y)dz + F [T* > t]H(y) 

where \R\ < 2F [W t E [x - A 1/2 ~ e/3 ,x]] < -J^A 1 / 2 "^ 3 < ^A £ / 3 since 

y/i > A 1/2 ~ 2e/3 . Finally, taking advantage of the estimates (13T31 and <KWi 
readily completes our proof. □ 

Proof of Lemmaim We take s = for notational simplicity. Introduce t a $ : = 
inf{t > : Vau t (A^)} and for 7 > write from Lemma and the 

notation of (13. 8ft (up to the same regularization procedure concerning F) 

F x [r A AT>A 20 } =PJ inf (f(0, x) + f V F(s, X?)a(<j>(a), X A s) )dW s 

+ ^(0,t l; x)) > 0,r A/3 > A 2/3+7 ] + O po ;(A) := Q, 



where under the assumptions of the Lemma, \Rp(0, t iy x)\ < Cti and F(0, x) < 
A a . For a given r > 0, consider the event A r = {3s < T : \X A — X^\ > r} 
where the increments of X A between two close times are large: by Lemma [HI 
it has an exponentially small probability. Hence, if we set 

M u := rVF(a s A )^(^),4))^ s := B {M)u , U = (M) u , 

J 

B is a standard Brownian motion (on a possibly enlarged probability space) 
owing to the Dambis, Dubbins-Schwarz Theorem, cf. Theorem V.1.7 in |RY99| . 
In addition, the above time change is strictly increasing on the set A C T and 
(M) t — (M) s > (t — s)a /2 (t > s) up to taking r small enough, because 
(A a -2) is in force. It readily follows that 

Q <¥ x [ o<i< mi^jM ti + CU) > A Q , T^ji > A 2/3+7 ] + O pol (A) 

<F X [ inf (B {i + 2Ca^U) > A a , r A , > A 2/3 + 7 , A^ + O pol {A) 
<F X \ inf (Br + 2Ca n 7 1 t i ) > -A a ,r A/3 > A 2/3+7 , 

inf (B s + 2Ca^s) < -A a ~<, A c r ] + O pol (A) 

0<s<(M) a2/3+7 

+ F x [t a p > A 2/3+7 , inf (B s + 2Ca^s) > -A""^], 
for C > 0. Thus, from Lemma M and standard controls 

Q < F x [3i : < i < A 2/3+7_1 , sup \B S - B- t . + 2C% \s - U) \ > A a " c - A Q 
r A , s > A 2 ^] + F x [ inf £ s > -A Q ^ - CA 2/3+7 ] + O poi (A) 

0<s<a A 2 ' 3 +^/2 

< O poi (A) + C(A a " c_/3 " 7//2 + A p+ ~i/ 2 ). 
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Choose now 7, ( s.t. (£ + |) = 77 > 0. The proof is complete. □ 

Proof of LemmaUR Taking also s = for notational convenience, we write 

P :=P,[r A < t, A^/ 2 F-(r A ,X T A A ) G [a, 6]] < O poi (A) 
L*/aj 

■ ^ « — 1 H— 1 

(3.15) 

using Lemma [8] for the last identity. 

A Taylor formula gives: F(t h X£) = F(U-i,X£_ J + ^(W^ - WVJ + 
Ht lA : = + #£_ 1)ti where S ti _, = VF a(U-i, X^_J , E^JIi^H < 
CA 2 . Conditionally to At ._ x has a Gaussian distribution 

^F(t i _ 1 ,X£jJE tw || 2 A). ' 

In addition, on the event X A _ l G Vao t ._ ( ro ), HS*^ || 2 A > aoA and we obtain 

-.^JF-^X*) G [aA 1/2 , 6A 1/2 ]] 
=^_ 1 [(M I _ 1 + J Rt 1 , t J"e[aA 1 / 2 ,6A 1 / 2 ]] 

<P^_ x G [-6A 1 / 2 - A 3 / 4 , -aA 1 / 2 + A 3 / 4 ]] 

+ P^_ 1 [|^_ 1 , t J>A 3/4 ,X A ^Aj 
^VJM^ G [-A 1/2 (6 + A 1/4 ), -A 1/2 (a - A 1/4 )]] 



+ CA 1/4 exp (-c- 



A 



using the Cauchy-Schwarz inequality and Lemma for the last inequality. 
Hence, we derive from (I3.15P 



P < g E,[l r A >ti _ liX A_ ieyaDti _ i(ro) (CA^exp j^- 



C A 



r AV"(-AV«, / (y-Fjt^XtJ)^ dy \ 

We now upper bound the above integral on the event {r A > t;_i} C {F(tj_i, X t A _ ) > 
0}. 

. If y < 0, clearly one has (y - F(t 4 _ l7 X A J) 2 > F 2 (t i _ 1 , X A J. 
. If y G (0, [A 1 /2(A 4 / 4 -a)] + ), one has (y-F^, X A J) 2 > fF 2 ^,^ J- 
^^iF^.XtJ-A^. 

Thus, we obtain that P is bounded by 

L*/AJ F 2 ft-_i X A } 

C(6 - a + A 1/4 ) E,[l rA>ii _ liX A 6y (ro) exp(-c ' ^ )] + O poi (A). 

., i — 1 z i — 1 / \ 

1=1 
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The end of the proof is now achieved by standard computations done in |GM04| 
p. 212 to 217. We only mention the main steps and refer for the details to 
the above reference. First, we replace the discrete sum on % by a continuous 
integral, then we apply the occupation time formula to the distance process 
(F(s, A^)) s < t a using the non characteristic boundary condition, as in the 
proof of Theorem 

(b-a + A 1/A ) r l F 2 (s X A ) 

P<C { - ± 1 E x [l T * >s>x . eVdDs{ro) exp(-c ^ s) )\ds + O pol (A) 

<C (& ~ a + Al/4) jT exp(-c^)E x [^ ArA (F(.,X A ))]^ + O pol (A). 

Then, we use (EU) to obtain P < C((b - a + A 1/4 ) which is our claim. □ 

Remark 13 Finally, we mention that if a a* is uniformly elliptic, the rest 
, + can be avoided and the result can be stated without the contribution 
A 1 / 4 . Indeed, we can directly exploit that the Euler scheme has conditionally a 
non degenerate Gaussian distribution and usual changes of chart associated to 
a parametrization of the boundary (see e.g. jGobOO^ ) give the expected result. 



4 Extension to the stationary case 

^.1 Framework 

In this section we assume that the coefficients in (11.11) are time independent 
and that the mappings 6, a are uniformly Lipschitz continuous, i.e. (X t ) t > is 
the unique strong solution of 

X t = x + f t b(X s )ds + f t a(X s )dW s ,t> 0, x e M. d . 
Jo Jo 

For a bounded domain D C R d , and given functions /, g, k : D — > R, we are 
interested in estimating 

u{x) := E x [g(X T )Z T + f 'f(X s )Z s ds], Z s = exp(- f S k(X r )dr), (4.1) 

Jo Jo 

where r := inf{t > : X t D}. 

Adapting freely the previous notations for Holder spaces to the elliptic setting, 
introduce for 9 e]0, 1]: 

(Ag) 1. Smoothness of the coefficients. b,a G Hi +( j. 
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2. Uniform ellipticity. For some a > 0, V(a?, £) G M d x M d , £*<7cr*(x)£ 

(D) Smoothness of the domain. The bounded domain D is of class H 2 . 
(Ce) Other coefficients. The boundary data g G H 1+e , f,kE H 1+ # and > 0. 

Note that under (Ag) and since D is bounded, Lemma 3.1 Chapter III of 
|Fre85| yields sup xg £, [r] < oo. Thus, (14.1 ft is well defined under our current 
assumptions. 

From Theorem 6.13, the final notes of Chapter 6 in |GT98| and Theorem 2.1 
Chapter II in Freidlin |Fre85| . the Feynman-Kac representation in our elliptic 
setting writes 

Proposition 14 (Elliptic Feynman-Kac's formula and estimates) 

Assume (A#), (D), (C#) are in force. Then, there is a unique solution in 
H 1+e n C 2 (D) to 

Lu — ku + f = 0, in D, 

(4.2) 

U\dD = 9, 

V 

(where L stands for the infinitesimal generator of X) and the solution is given 

by dHD. 

In the following we denote by F(x) the signed spatial distance to the bound- 
ary 3D. Under (D), D satisfies the exterior and interior uniform sphere con- 
dition with radius r > and F G H 2 (VaD(r )) where Vqd{tq) := {x G M. d : 
d(x, 3D) < r }. Also, F can be extended to a H 2 function preserving the 
sign. For more details on the distance function, we refer to Appendix 14.6 in 
|CT98| . 

4-2 Tools and results 

Below, we keep the previous notations concerning the Euler scheme. We also 
use the symbol C for nonnegative constants that may depend on D, b, a, g, f, k 
but not on A or x. We reserve the notation c for constants also independent 
of D,g,f,k. 

We recall a known result from Gobet and Maire |GM05| (Theorem 4.2) which 
provides an uniform bound for the p-th moment of r A : 

Vp > 1, limsupsupE x .[(r A ) p ] < oo. (4.3) 
Let us now state the main results of Section [2] in our current framework. 
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Proposition 15 (Tightness of the overshoot) Assume (Ag-2), and that 
D is of class H2. Then, for some c > 0, 

supE :E [exp(c[A" 1/2 F^(X^A)] 2 )] < +00. 

A>0 

From the proof of Theorem [3] and the estimate (14.31 ) we derive: 

Theorem 16 (Joint limit laws associated to the overshoot) Assume (Aq), 
and that D is of class H2 . Let tp be a continuous function with compact support. 
With the notation of Theorem^ for all x G D, y > 0, 

ZMX T )(l - H(y/\VFa(X T )\))~ . 

4-3 Error expansion and boundary correction 

For notational convenience introduce for x 6 D, 

u(D) = E x (g(X T )Z T + f Z s f(X s )ds), 

Jo 

u A (D) = E x (g(X A A )Z A A + Z A (s) f(X A s) )ds). 

J 

The second quantity is well defined owing to (14.31) . 

Theorem 17 (First order expansion) Under (Ag), (D), (Cg), for A small 
enough and with the notation of Theorem^ 

Err(A,g,f,k,x)=u A (D)-u(D) 

= coVAK x (Z T (Vu - Vg)(X T ) ■ VF(X T )\VFa(X T )\) + o(VA). 

Define now D A = {x e D : d(x,dD) > c y/A\VFa(x)\}. Introduce f A = 
M{ti > : X A e D A }. Set 

u A (D A ) = E x [g(X A .)Z A A + jT* Z A s) f(X A s) )ds}. 

One has: 

Theorem 18 (Boundary correction) Under (Ag), (D). (Cg) and assum- 
ing additionally VF(.)\VFo~(.)\ is in C 2 , then for A small enough one has 

u A (D A ) — u(D) = o(v / A). 



E x [Z a MX A a)1 f ~ {x * te^] — > E 
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4-4 Proofs 



Note carefully that all the constants appearing in the error analysis for the 
parabolic case have at most linear growth w.r.t the fixed final time T. Estimate 
(14.31) allows to control uniformly the integrability of these constants in our 
current framework. Thus, since the arguments remain the same, we only give 
below sketches of the proofs. 

Proof of Propositioning It is sufficient to prove that there exist constants c > 
and C s.t. VA > 0, sup A>0 F x [F~(Xf A ) > AA 1 ' 2 ] < Cexp(-5A 2 ). Then any 
choice of c < c is valid. For x G D, we write 

P:=P X [F-(X T A A ) > AA 1 / 2 } 

= EEfW^l^ <u nF-(X^)>AA^ T . ]] 
. ^ T i %— i l i— i 

i€N* 

where t£_ := inf{s > : X A ^ D}. From Lemma [HI we get 

P < Cexp(-dA 2 ) J2 > fc-i,^ < 
ieN* 

Lemma 16 from |GM04| remains valid under our current assumptions and 
yields 

P < Cexp(-5A 2 ) £ Eil^^MX* <£ D\F ti _A + O pol (A))}. 
ieN* 

On the one hand, £ieN* ^■T A >t i - 1 ^-x A (D = 1t a <oo = 1 owing to (14.31) . On the 
other hand, we have £ feN *Pz[T A > U-i] = A^E x [t a ] < C/A using 
again. Finally, we obtain that P < Cexp(— cA 2 ) which concludes the proof. □ 

Proof of Theorem [73 Similarly to the proof of Theorem [6] we suppose first 
that u G H 3+ 0. The general case can be deduced as in the parabolic case using 
suitable Schauder estimates, given in the final notes of Chapter 6 in |GT98| . 
see also our Appendix. 

In this simplified setting, keeping the notations introduced in the proof of 
Theorem [6j we obtain 

Vii{A )g JXx) = Z%{Vu-Vg){ii dD {X^))VF{X^)F-{X^) 

+ (^1 1i<t a[i /1 ,0(A 1 ?) (4.4) 
ieN 

+l un cl r . . i (1 eK (w(X A )Z, a -m(X a )Z, a 

+Z A /(X A )A)l)v 0>rA . (4.5) 
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Since the constant in (12.61) depends linearly on time, the contribution associ- 
ated to the remainder (14.41) can be bounded by CA ^ x (A~ 1 E :r [r A ]). From 
(14.31) . this quantity is a 0(A + 2 e ) = o(A2) for e small enough. Similarly to 
( 12.101) the term (14. 5fl can be bounded by 

E[(^l ti<T Al ( ^ )C 0(A 2 {l + luloo + IVmU + |D 2 u|oo + |D 3 u|oo} 

+ A^[D 3 uU))]+O pol (A) 
<CA i ^E[r A ]=o(A 1 / 2 ). 

We eventually derive the result as in Section [21 □ 

Theorem [18] can be proved as Theorem using a sensitivity result analogous 
to Theorem 2.2 in |CGK06| for elliptic problems, see e.g. Simon [SimSOj. We 
skip the details. 



5 Numerical results 



The numerical behavior of the correction of Theorem [5] had already been 
illustrated for the killed case in Section 3 of |Men06| . Additional tests are 
presented in |Gob09| . We now focus on the stopped case with the following 
example. Take d = 3 and introduce the following diffusion process 



dX t = b(X t )dt + a(X t )dW t: Vx e R 3 , b(x 
( 



a(x) 



M) 1/2 

-N) 1/2 (!) 









(x 2 x 3 Xt)* , 



l(i + k 2 i) i/2 (i) 1/2 (i + \x2\y /2 ) 



1/2 



a/2 



(5-1] 



and X to be specified later on. Set D = B(0,2). We consider an elliptic 
problem. Starting from a given function u(x) = X1X2X3 defined on D, we derive 
the PDE of type (14.21) associated to ( 15.11) satisfied by u by taking g = u\qd, 
setting / = —Lu where L stands for the infinitesimal generator of X in (15.11 ) 
and k — 0. One can easily check that —f(x) = x\x 3 + x\x\ + x\x 2 + f F3(l + 

\x 1 \)y 2 {l + \x^\) 1 ' 2 +x 1 (|) 1/2 (l + |rci|) 1 / 2 (l + |a: 2 |) 1 / 2 ]. Thus we have an explicit 
expression for the solution of (14.21) . 

For Xq s.t. (xq)i<j<3 G {—.7, —.3, .3, .7}, we take N MC = 10 6 sample paths 
for the Monte Carlo simulation and let A vary in {.01, .05, .1}. For all the 
computations, the size of the 95% confidence interval always varies in [1.5 x 
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10 -3 , 2 x 1CT 3 ]. For the absolute value of the absolute and relative errors over 
the 3 x 4 3 = 192 points of the spatial grid, we report the results in Table [TJ 
These results for the correction seem to indicate that the remainder o(A 1 / 2 ) 
in Theorem d8] is actually a 0(A). This will concern further research. 



A 


Without correction 


In the corrected domain 


.1 


0.169 (199%) 


0.0220 (24.4%) 


.05 


0.114 (133%) 


0.0115 (13.1%) 


.01 


0.0471 (54.7%) 


0.0026 (2.98%) 



Table 1 

Supremum of the absolute error for the Euler scheme (relative error in % in paren- 
thesis) 

In Tables [2] and El we also report the results obtained for the spatial points 
x = (-.7, .3, .7) and x = (-.7, .7, -.7). 



A 


Without correction 


In the corrected domain 


.1 


-.0913+/- .0019 


-.1477 +/- .0016 


.05 


-.1051 +/- .0018 


-.1465+/- .0016 


.01 


-.1282 +/- .0017 


-.1476+/- .0016 



Table 2 

Estimated value at xo = (—.7, .3, .7) (with 95% confidence interval). True value 
u(x Q ) = -.147. 



A 


Without correction 


In the corrected domain 


.1 


.5368 +/- .0019 


.3866 +/- .0016 


.05 


.4648 +/- .0018 


.3634 +/- .0016 


.01 


.3851 +/- .0016 


.3473 +/- .0016 



Table 3 

Estimated value at xq = (—.7, .7, —.7) (with 95% confidence interval). True value 
u(xq) = .343. 



Eventually, for the Monte Carlo method, taking x = ( — -7, .3, .7) and the 
previous values of A, in Figure[5]we plot — log(ErrA/c) m function of — log(A), 

where Err MC := (g( x %) + J q f( x f(s)) ds ^ ~ u ( x o)- The curve 

is quite close to a right line with slope 1/2 as it should from Theorem [T7l 
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2 2.5 3 3.5 4 4.5 5 

-Log(Delta) 



Figure 5. Error for the Monte Carlo method (without correction) as a function of 
A, in logarithmic scales. Evaluation at xq = (—-7, .3, .7). 



6 Conclusion 



We have proposed and analysed a boundary correction procedure to simulate 
stopped/ killed diffusion processes. This is valid for non-stationary and station- 
ary problems, in time-dependent or time-independent domains. The resulting 
scheme is elementary to implement and its numerical accuracy is very good 
in our experiments. The proof relies on new asymptotic results regarding the 
renormalized overshoots. 

To conclude, we note that the boundary correction procedure is very generic 
and could be at least formally extended to general Ito processes of the form 
dX t = b t dt + a t dW t . In that case, the smaller domain would be defined lo by 
uj replacing VF(t,x)a(t,x) by VF(t,X t )a t . Even if our current proof relies 
on Markovian properties, we conjecture that the correction should once again 
give a o(\fK) independently of the Markovian structure. Numerical tests in 
|Gob09| support this conjecture, which will be addressed mathematically in 
further research. 



A Proof of Theorem [6] in the general setting 



In this section, we detail how the proof of Section [2] has to be modified under 
the assumptions of Theorem [U i.e. for g G H 1+e and without compatibility 
condition so that u G H 1+e . Actually, u is smooth inside the domain but high 
order derivatives may explode close to the boundary. These features have to 
be accurately quantified to show that the induced singularities are integrable. 
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A.l Preliminary notation and controls 



Introduce the parabolic distance pd: for (s,x), (t,y) G V, pd((s,x), (t,y)) = 
max(|s — t\ 1 / 2 , \x — y\). We also denote for a closed set A G V and (s,x) G 
V, pd((s, x),A) the parabolic distance of (s, x) to A. Note that pd((s, x),Wf] 
{v > s}) > min(F(s,x), — s), so that we obtain the easy inequality: 

<TT7— T + ^^- (A.l) 



pd((s,x),PI?n > s}) - F(s,x) v/T^7 

Under our current assumptions, for some constant C > 0, we have 

|D 2 m(s,x)| + \D 3 u(s,x)\ < Cpd{(s,x),VVn{v > s})~ 2 ; (A.2) 

for {t,y) / (s,z), — — 

pd((s,x), 

< C[pd((s, x),VVn{v> s}) A pd((t, y),W n{v > t})]~ 2 ~ e ; (A.3) 

fort ^ s ' |t- 3 |(i+»)/g 

< C[pd((s, i),PPn > s}) A pd((t, x) 1 VVf]{v> t})}- 2 - 9 . (A.4) 



The above constant C is uniform w.r.t. (s,x) G V,(t,y) G P or G 
P. These inequalities are obtained with the interior Schauder estimates for 
the PDEs satisfied by the partial derivatives (d x .u)i<i<d, see Theorem 4.9 in 
[Lie96] . 

We first state an important proposition for the error analysis with possibly 
explosive controls as in ( 1A.2I )- (IA.3I )- (IA.4I) for the derivatives. Namely, under 
our current regularity assumptions, in order to perform a Taylor expansion 
we have to work with interior points located in small balls, which distance 
to the boundary is uniformly bounded from below within the ball. The next 
proposition states that this is the case if the ball centers are "far enough" from 
the side of V. 

Proposition 19 Assume V G H 2 and take e G]0,1[. For all (t,x) G V fl 
V dv {r /2)\V av {2^ 1 l 2 ^) (r is defined in Section EXf , one has for \/y G 
B(x, A 1/2(1 " £) ) and s G [t,t + A] 

F(s,y)>^F{t,x) 

for A small enough (uniformly in t,x, s,y). In particular, y belongs to D s . 
Proof. Since F G H 2 , one has 

F(s, y) >F(t, x)-CA + (VF(t, x),y-x)- CA 1 " 6 . 



35 



The norm of WF(t, x) equals 1, since VF(t,x) is the unit inward normal 
vector at the closest point of x on dD t . Therefore, for A small enough and 
using \F(t,x) > A^ 1 "^, we have 



F{s,y) > F{t,x) - -A^~ £) > -F{t 



2" -4-- Xj ' 
which is the expected inequality. □ 

We are now in a position to deduce useful local upper bounds for the deriva- 
tives of u and their Holder-norms, under the assumptions of Theorem [6l 

Corollary 20 Take e e]0, 1[. There exists a constant C > such that for A 
small enough, for all (t, x) E V\V dv i2A l l 2 ^) , for all (y, z) E B(x, A 1 ' 2 ^) 
and (r, s) E [t,t + A], we have 

\D 2 u(s,y)\ + \D 3 u(s,y)\ < j^—j + ^- t ; (A.5) 

\D 3 u(s,y)-D 3 u(s,z)\ ^ C C 
f ° ry ^ Z > [y—A e -F^(^) + (T-t)W ; (A - 6) 

\D 2 u(r,y)-D 2 u(s,y)\ C C 

J0T T ' | r _ s |(i+0)/2 "F 2 + e (t,x) (T-t)W 1 ' 



Proof. Note that if (t,x) E D\14 c (2A 1 / 2 ( 1 - £ )), we have T — t> 4A 1 ~ £ . 
Estimate $A.5\) . In view of (IA.2I) and (lA.ljl . the upper bound of \D 2 u(s,y)\ + 



D 3 u(s,y)\ is equal to F2 F > + On the one hand, by easy computations, 



F 2 (s,j/) 1 T-s 

we prove 



1 1 T-t 1 1 C 
< < < 



T-s ~ T-tT-t-A ~ T-tl- A e /A ~ T-t 
for A small enough. On the other hand, we have 

1 < ° 



F(s,y) ~ F(t,xY 

Indeed, if x is far from D t (and thus y far from D s ), both terms F(s,y) and 
F(t, x) are bounded from above and from below. In the other case when (t, x) E 
VnV av (r /2)\V dv (2A 1 / 2 ^), Proposition M yields > \. Therefore, the 
upper bound (1A.5I) readily follows. 

Estimates ( L4.6|] and ( L4. 7|J . They are proved following the same arguments, 



the details of which are left to the reader. □ 
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A. 2 Error analysis 



Recall from the previous proof of Theorem [6] that the main term to analyse is 



e 22 — 



( E * 

0<U<T 



u(t i+ll x*)z£ 



-U(U, X*)Z$ + Zgffa X£)A\) l T r 0>T A 



^ • • •)l r ''0>T A AT + ( ' ")lr r o>T A AT ; — + e 222) 

0<t,<T-4A 1 - £ T-4A 1 - £ <t 4 <T 

where we have just splitted the summation on tj. 

Control ofe^i- The idea is to perform a stochastic expansion of u(t i+ i, X£ )Z£ 
—u(ti, X£)Z£ + Z^f(ti, X^)A as in (12.8ft . Under our current assumptions, the 
difference comes from the high-order derivatives that are no more uniformly 
bounded or uniformly Holder but only locally, with local estimates given in 
Corollary EHl Thus, following the same computations that have led to (12.101) . 
we obtain 



e 221 — 



lt^l^cl x*eB(x?A? (1 - £) ) 

T-iAl-E E » 

1 1 



0<t i <T-4A 1 - £ 

. 2 , A i ^A 



o( ( A^A|xf-x t tn(^ )Xtt) , T _ t 



O (A 



+ Mti,*i+l )lr''0>r A AT- 



(A.f 



The derivatives appearing in (Mt j ,t j+1 )o<t i <T (see equations (12.91) and (12.101) ) 
are controlled by (1A.5I) on (Af.) c . The control of (|2.1ip remains valid for the 
(M tuti+1 ) Q < ti<T that yields a negligeable contribution. It follows that 



14,1 <CA^ 



E 



0<t i <T-4A 1 - e 



1 U<r Al F{t ^ X A ) y 2A ^l-,) A 



1 



Standard computations show that 

A 1+8 A i+e /•T-4A 1 -+A rft 

A" > < A~ / 

0<i j <T-4A 1 -< s ^ fc V ,yU 

which implies 



1t ? '0>t a AT- 



( T _ t ) 1+0/2 



0(A 



life. 

2 ' 2 



e* 21 | <CA^ 



TAt 



F(^( t)) xA ) e[ 2AV2Ci-e) )r . o/2] F(0(t),X^ (t) )- 2 -^t]+O(At + ^ 
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Adapting the previous analysis of Section [2] for the term e^ n , we get 

14,1 <CA^n TAT l F{t . xmA >, v -,>.3r l ,, i )F(t,Xf)- 2 - e dt) +0(A*+*) 
< CA ! ?p_ i)9 -^(F(.,X i )) i i ! ,) + 0(Ai+% 

using Lemma[8]for the last but one inequality, and the occupation time formula 
for F(t,X^) for the last one (recall that a is uniformly elliptic). 
Finally using ( 12.61) . one gets 

I4 x | < ca^ {^f^ / l ^ 1 _ E) y~ 2 ~ e (v + A 1/2 )<iyj + o(A^ + ^) < caH = (av*). 

Control ofe^- Apply a Taylor formula with integral rest at order one in space. 
The ^-Holder continuity in space of Vii and the (1 + 6>)/2-H61der continuity in 
time of u directly give a contribution in 0(A 1 ^ 2+9 ^ 2 ~ e ) = o(A 1 / 2 ) for e small 
enough. This completes the proof. □ 
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